clc;
clear;
syms f z u v x y v1 v2 f0 f1 f2;
u = x ^ 2 - y ^ 2 - x * y;
v1 = int(diff(u, x), y);
v2 = int((diff(u, y) - diff(v1, x)), x);
v = v1 + v2;
f0 = u + 1i * v;
f1 = subs(subs (f0, x, z), y, 0);
f2 = (-1 + 1i) - subs(f1, z, 1i);
f = f1 + f2